Enhancement of the electric dipole moment of the electron in PbO 
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The a(l) state of PbO can be used to measure the electric dipole moment of the electron d e . We 
discuss a semiempirical model for this state, which yields an estimate of the effective electric field 
on the valence electrons in PbO. Our final result is an upper limit on the measurable energy shift, 
which is significantly larger than was anticipated earlier: 2|Wd[d e > 2.4 x 10 25 Hz [^^] . 
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In his pioneering work, Sandars pointed out that the 
effective electric field on a valence electron in a heavy 
atom is enhanced by a factor ^a 2 Z 3 relative to the ap- 
plied laboratory field B. That started a long search for 
the electric dipole moment (EDM) of the electron d e in 
atomic experiments (0). The most stringent limit on d e 
follows from an experiment on atomic Tl (Z=81) |J. 

Even larger enhancement is present in heavy polar di- 
atomic molecules @, || . The heavy atom there is sub- 
jected to an internal E-field of ~1 a.u. w 5 • 10 9 V/cm, 
which is further enhanced by the relativistic factor a 2 Z 3 . 
This effective field is many orders of magnitude larger 
than available laboratory fields; this makes diatomic 
molecules very attractive systems to look for d e . 

Since d e is linked to the electron spin, one must work 
either with radicals, which have an unpaired electron 
in the ground state, or with excited states of "normal" 
molecules. Diatomic radicals with the ground state £1/2 
have large enhancement factors which can be relatively 
easily calculated ||, J5J . The first results of an EDM mea- 
surement in such a molecule (YbF) were recently pub- 
lished H |i| . The molecule PbO is a favorable candidate 
for a search for d e in the excited state a(l) ^, |o) , and the 
group at Yale has begun EDM experiments on PbO . 
It is therefore timely to estimate the effective internal 
field for the state a(l) of PbO. 

The interaction of d e with an electric field E can be 
written in four-component Dirac notation as ]T^ ]: 




crE 



(1) 



After averaging over the electronic wavefunction this in- 
teraction can be expressed in terms of an effective spin- 
rotational Hamiltonian || ^| : 

re// 



H 



W d d e (J e ■ n), 



(2) 



where J e is the electronic angular momentum and n is 
the unit vector along the molecular axis. In this paper 
we estimate Wd for the molecular state a(l): 



W d = d- l {a{l)\H d \aO-))i 



(3) 



where we used that (a(l)| J e -n\a(l)) = fi[a(l)] = 1. The 
f2-doubling for states with 17 = 1 is very small and even 
in a weak external electric field the energy eigenstates 
correspond to definite fl rather than definite parity. The 
energy of the molecule can be then written as: 



W(J, M, fi) = BJ(J + 1) + W d d e n 



DEqUM 
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(4) 



where B is the rotational constant, Eq is the external 
electric field, and D is the molecular dipole moment. The 
EDM contribution can be determined from the difference: 



W(J, M, n) - W(J-M-ft) = 2W d d e 



(5) 



In order to estimate the matrix element in Eq. (||) we 
construct here a semiempirical wavefunction of the state 
a(l). We use the MO LCAO approach, where each molec- 
ular orbital is expressed as a linear combination of atomic 
orbitals, and all molecular matrix elements are reduced 
to the sums of atomic matrix elements. The HFS or SO 
interactions as well as the EDM enhancement factor grow 
very rapidly with nuclear charge Z . Therefore, we are in- 
terested only in the Pb part of the MO LCAO expansion. 

Analysis of the molecular observables requires knowl- 
edge of several atomic matrix elements for Pb. We cal- 
culate these in the Dirac-Fock approximation both for 
neutral Pb and for Pb + . Results are given in Table [| for 
the orbitals 6s and 6pj. For the HFS operator we calcu- 
late the parameters hk,k' as defined in H (we use atomic 
units unless the opposite is stated explicitly): 



hk,k 



2m v 



{fk9k> +9kfk')dr, 



(6) 



where g n = 0.59 is the nuclear g-factor of 207 Pb, m p is 
the proton mass, fk and gk are upper and lower compo- 
nents of the Dirac orbitals, and k = (I — j)(2j + 1) is the 
relativistic quantum number. For the EDM operator (|l|) 
in our minimal basis set there is only one nonzero radial 
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TABLE I: Atomic parameters h k ^y (GHz), w sp (a.u.), and £ 
(cm -1 ), calculated in the Dirac-Fock approximation for Pb 
and Pb + . The relativistic quantum number k is equal to —1, 
1, and —2 for Si/ 2 , P1/2, and p 3 / 2 correspondingly. 







h-2,-2 


W sp 








Pb + 45.5 


-8.9 -1.1 


1.8 


-34.1 


-1.7 


9452 


Pb 42.3 


-7.5 -0.9 


1.4 


-30.0 


-1.7 


8077 


ratio 1.08 


1.20 1.25 


1.30 


1.14 


1.0 


1.17 



integral, between 6S1/2 and 6pi/2 orbitals: 



d(f> 2 
g-igi^-r dr, 
or 



(7) 



where </> is the atomic electrostatic potential. Finally, we 
also need the atomic SO constant £ for the 6p shell: 



#so = ^ • s > 



£ - g ( £ 6p 3 /2 _ £ 6p 1/2 )- (8) 



The last row of Table § gives the ratio of the radial 
integrals for the ion and for the atom. Note that the 
ratios are similar for all relevant integrals, and a simple 
relation between w sp and HFS constants holds for both 
cases: 



1.7^-/1-1-1/11,1. 



(9) 



This relation is critical for semiempirical models of the 
EDM enhancement: it implies that the value of the EDM 
enhancement does not depend on what set of radial inte- 
grals (i.e. atomic or ionic) is used. Similarly (|TJ) holds for 
other principal quantum numbers n, e.g. for 7s and 7pj. 
The particular choice of the radial integrals becomes im- 
portant for the normalization of the wavefunction. The 
Pb atom in PbO is positively charged; therefore, below 
we use the ionic set of integrals from Table ^[ 

In order to develop a semi-empirical model for the state 
a(l), we have found it necessary to also consider the 
wavefunctions of several low- lying states of PbO. Pre- 
vious work has shown that these states correspond to 
the configurations and nominal A, S-coupling terms as 
follows rfL3 , 



a(l)[ 3 S+],C(0 + )[ 3 So 

J 4(o+)[ 3 n ],s(i)[ 3 n ] 



X(0+)[%+] 
,C'(1)[ 3 £7] 



,.,2 2 2 4 
1 2 3 1 



(10) 

^a|a|7r 3 7T2,(ll) 
ofcr|o- 3 7r^7r 2 ,(12) 



In the naive ionic model of PbO, the two 6p electrons 
from Pb move to O and close its 2p-shell. This suggests 
that the orbitals 171,2 and 7Ti are centered on O, and that 
03 is predominantly of the Pb 6s-type. However, below 
we do not impose any constraints on the MO LCAO coef- 
ficients for these molecular orbitals, based on this expec- 
tation. Note that only the orbitals 03 and 7Ti 2 contribute 



to the spin-density of the molecular states under consid- 
eration. Thus, we do not need to know the orbitals 01,2, 
and below we omit the index for the orbital 03. 

Now we specify coefficients of the MO LCAO expan- 
sion for the three valence orbitals of interest: 



& 1 6s 1/2 
I 



2^/ ^|6pi/ 2 ,o.) -+ 

p i (2wy||6p 1/ 2,, 



|6p 3 /2 



|6p. 



3/2, [ 



Pi\6P3/2 



(13a) 

,(13b) 
(13c) 



where co = ±1/2 and oj' — ±3/2. The numerical coef- 
ficients are chosen to account for the quantum number 
A: A = 0(1) for a(n) orbitals. In order to calculate Wd 
we must determine the 4 parameters in Eqs. (|l3|). Below 
we try to constrain these parameters using experimental 
information about states (11) and (fl2|). To simplify the 
notation, we define the ground state of the molecule as 
vacuum. Then each of the excited states in Eqs. ( pT| ) 
and ( |l2] ) is a two-particle state with one hole and one 
electron. We do not use any special notation for the hole 
states; instead we simply write the hole orbital in front 
of the electron one. We construct wavefunctions of these 
states from the orbitals (jp|), using at the first stage the 
A, S-coupling scheme classification: 

Ki)[ 3 £+]} = 

= -^(|7Tl,A=-17r2,A=l) + |7ri,A=l7r 2 ,A=-l)) | TT) 

= 7^= (ki, -1/2712,3/2) + ki, 3/2712,-1/2)) ; (14a) 

|A(o+)[ 3 n ]) = 

= ^(I^,a=i)|U)-|^,a=-i)|TT)) 

= -7= (1^-1/271-2,1/2) - |cti/ 2 7T2 -1/2)) ; etc. (14b) 

We write each wavefunction in both A— a and lu—lu repre- 
sentations; the latter is more convenient for our purposes. 

The rules for calculating hole matrix elements follow 
from the fact that the hole in the state \u>) actually means 
the absence of the electron in the state \ —u>). Thus, the 
expectation value for an electronic operator P over the 
hole state \u) can be written as: 



(u\P\w) 



-w\P\-w) b =t(u\P\u) 



(15) 



where we applied the time-reversal operation T . Thus the 
final sign depends on the time-reversal symmetry of P, 
with the minus sign corresponding to a T-even electronic 
operator. For example, the HFS interaction is given by 
the product of the T-odd electronic vector AJ e and the 
nuclear spin I. Thus, for the HFS interaction the plus 
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sign in Eq. ([15]) is correct. A similar argument shows 
that the SO constant £ for a hole has the opposite sign 
as for an electron. 

From Eqs. (jl^), the first-order SO splitting A AB be- 
tween states ^4(0 + ) and S(l) is: 



&AB = ^ [( 7r 2,3/2|'s|7r2,3/2)-(7r2,l/2^s|7r 2 ,l/2)] = 



Using the experimental value of this splitting |15| and 
the ionic value for £ from Table H, we estimate P 2 : 



P? = 



2A AB _ 2 • 2420 
£ ~ 9450 



= 0.51. 



(17) 



We see that the orbital 7T 2 has a large contribution from 
the Pb orbital 6p. The data on energy levels [HJ shows 
then that for all levels with one electron in the 7T2 or- 
bital, the SO interaction is comparable to the splittings 
between these levels. Therefore, there must be significant 
SO mixing between such states. 

We start with the mixing within configuration o 2 'k\'Ki- 
The mixing angle a between states <x(l) and C"(l) is: 



5 E 1 -| J ff so | 3 s+) a +6 



»aC" 



2IA. 



aC 



(18) 



where A. a c> is the energy splitting between a(l) and 
C"(l), and & = £P?. If we assume that Pf < P 2 2 (corre- 
sponding to the naive ionic model), we can estimate the 
value of a: 



a 



6 



0.3, 



2|A aC v| 

and write the new wavcfunction in the form: 



(19) 



|a(l)) = C a Kl,3/27T 2 -1/2) + Sakl -1/2^2,3/2), (20) 

c a = cos ^— — , s a e sin ^— — a^j . (21) 

SO interaction also mixes configuration cr 7rf 7T 2 with 
configurations ait^ir^ and oi:\-k 2 . These mixings can 
be accounted for by substitution of the original orbitals 
Ki,i/ 2 ) with the perturbed orbitals 



Fi,l/2/ 



CiFi,l/2> + »ikl/2> 



(22) 



There is no experimental information about levels of the 
configuration ait^ir 2 , so we cannot reliably estimate the 
mixing parameter s 2 . In contrast, both levels with f2 = 
1 of the configuration aTvfiv 2 are known [i.e., B(l) and 
D(l)]. That allows us to write for si the estimate: 



2.8s 2 a PiS p . 



(23) 



These SO mixings then lead to the final form of the wave- 
function of the state a(l): 

|o(l)> = £^1^,3/2^2,-1/2) + s a 1*1-1/2^2,3/2)- (24) 



The G-factor for the state (E4T) is given by: 



G|, = (a(l)|L + 25o|o(l)) = 2-s 2 a s( 



2 „2 2 „2 



(25) 



The measured value Gu — 1.84(3) Q corresponds to the 
following equation for mixing parameters: 



4*1 + 4*2 = 0.16(3). 



(26) 



(16) The signs of the parameters si^ should be chosen so that 



the contribution of atomic orbital 6pi/ 2 to the molecular 
orbital a is increased: in this case relativistic corrections 
to the binding energy of the a orbital are positive. 

The matrix element of the HFS interaction for the state 
a(l) @ has the form: 

{a(l)\H Ms \a{\)) = (27) 

= 4 [( 7r l,3/2khfsKl,3/2) - {^2,1/2 1 ^hfe 1^2,1/2)] 
+ 4 [( 7r 2,3/2khfsk2 : 3/2) ~ (*l,l/2 1 hhis 1*1,1/2)] • 

We use expressions from Ref. || for the one-electron ma- 
trix elements and numbers from Table [|, combined with 
the measurement of the hyperfine constant for the state 
a(l): ^4 1| = —4.1 GHz ^H^, to find another equation re- 
lating the various coefficients of the model: 



30( C 2 S 2 + S 2 S 2^2 + 18 ( C 2 S 2 + 

+ (4.6<£c? - 1.44) Pi +(4-64c2 - 1-44) P 
-A.7slciSiPiS p -4.7c 2 a c 2 s 2 PiS p = 4.1. 



4 s ?) Sp 

2 ^ 



(28) 

(Note that the formulae of Ref. || are strictly applicable 
only for orbitals and states with u> = Q = 1/2. Eq. ( psj ) 
takes into account simple modifications of these formulae 
for the present situation.) 

We now have five equations, namely: dT^), (|lS|), (p3), 
(26), and ( |28| ) on seven parameters of the wavefunction 
(24). That leaves us with two independent parameters 
of the model. We introduce two additional constraints, 
which account for normalization and the Pauli principle: 



si 



Pt 



Pi < N . 



(29) 



We choose No = 1.2 here in order to account for inac- 
curacy of the Hartree-Fock approximation used to deter- 
mine the atomic parameters in Table |. 

The parameters a and P 2 are unambiguously fixed by 
Eqs. @ and |TJ). We choose s± and Pi as free param- 
eters and solve Eqs. (|23|), (26), and ( p8|) for parameters 
S2, S p and S s . After that we reject solutions which do 
not meet the constraints (f29|). A typical solution is: 



a = 0.3; 
si = 0.107; 
s 2 = 0.449; 



P 2 = 
Pi = 
Sp = 



0.714 
0.503 
0.349 



(30) 



S„ = 0.812. 



Some of the parameters appear relatively well-defined, 
while others are not. The variation ranges are: 



si 
S:, 



0.2; 
0.5; 



0.4 < 



S 2 

S 2 



0.5 
0.5. 



(31) 
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TABLE II: Dependence of the EDM constant W d (in a.u. 
the parameters of the model. 





G|| 


a 


Pi 


\w d \ 


(GHz) 












—4.1 


1.84 


0.30 


0.51 


19.6 


13.7 


-4.1 


1.81 


0.30 


0.51 


19.1 


12.0 


-4.1 


1.87 


0.30 


0.51 


20.3 


15.4 


-4.1 


1.84 


0.24 


0.51 


19.8 


13.2 


-4.1 


1.84 


0.36 


0.51 


20.3 


15.5 


-4.1 


1.84 


0.30 


0.41 


18.2 


12.1 


-4.1 


1.84 


0.30 


0.61 


20.3 


15.1 


-3.3 


1.84 


0.30 


0.51 


17.0 


11.2 


-4.9 


1.84 


0.30 


0.51 


22.0 


16.4 



The parameter P\ appears to be restricted only by the 
normalization condition (^9|). 

It may be possible to add some restrictions to reduce 
the ranges of variation in Eq. (|3l]). For example, the 
relatively large value of S2 should require a large value 
of S p . However, such additional restrictions would add 
arbitrariness to the model and may affect its reliability. 
We use only the minimal set of constraints to determine 
the range of possible values of Wd- 

For the wavefunction ( |24| ) , there are two contributions 
to the EDM parameter Wd from each of the one-electron 
orbitals with \u\ = 1/2: 



W d = 



-c 2 W? 2 

u a*' d 



:WJ! 



4w 



Sp 



V3 



(V2c 2 a C2S2P2 



2 2 o 
c a s 2 D P 



1 s 2< ? 



(32a) 



(32b) 



We find that the first term in ( |32b| ) always dominates 
the sum. The second term is not negligible, but the final 
two terms contribute < 10%. It is important that the 
leading contribution to Wd is similar to the first term in 
Eq. (Pq), which dominates the HFS. This implies that the 
parameter Wd is well-constrained even though some of 
the parameters of the wavefunction are not. We obtain: 



IWrfl = 16.6 ±3.0 a.u., 



(33) 



where the uncertainty reflects the range of values found 
within the model just described. 

It is also important to check how Wd depends on the 
"fixed" parameters a and P 2 , as well as on the input data 
for A|| and Gy, since our model relating the MO LCAO 
coefficients to these parameters is rather crude. In Ta- 
ble H we solve the model equations for values of these 
quantities varying from the best values by ±20%. We 
find that this variation of the input parameters widens 
the range for Wd substantially (to ±5.4 a.u.), but still 
does not allow dramatically smaller values of W d - 

It is known from previous calculations of Wd for other 
diatomic molecules, that correlation corrections tend to 



decrease the result by 10-20% from the Hartree-Fock 
level. Therefore, we state our final result as a conser- 
vative lower limit on Wd- 



\W d \ > 10 a.u. = 12- 10 2 



Hz 



e cm 



(34) 



This lower bound is several times larger than earlier, 
naive estimates which did not consider the effect of SO 
mixing on the (nominally) 7r-type orbitals of the a(l) 
state H . Our model shows significant similarity between 
the orbital 712,1/2 i n PbO and the single valence orbital in 
the ground state of the free radical PbF. It is thus nat- 
ural that our bound is close to the value calculated for 
PbF ^T|. (Coincidentally, our bound is also similar to 
the calculated value for YbF 0, [If, g |||.) However, 
we stress that this first semiempirical estimate of the ef- 
fective field in PbO has very limited accuracy. Thus, 
more elaborate calculations of the a(l) state are highly 
desirable. 
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